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ABSTRACT 

We present new radial velocities from Keck Observatory and both Newtonian and Keplerian solu- 
tions for the triple-planet system orbiting HD 37124. Th e orbital solution for the system has improved 
dramatically since the third planet was first reported in lVogt et al.l (|2005f ) with an ambiguous orbital 
period. We have resolved this ambiguity, and show that the outer two planets have an apparent period 
commensurability of 2:1. A dynamical analysis finds both resonant and non-resonant configurations 
consistent with the radial velocity data, and constrains the mutual inclinations of the planets to be 
<- 30°. 

We discuss HD 37124 in the context of the other 19 exoplanetary systems with apparent period com- 
menserabilities, which we summarize in a table. We show that roughly one in three well-characterized 
multiplanet systems has a apparent low-order period commensuribility, which is more than would 
naively be expected if the periods of exoplanets in known multiplanet systems were drawn randomly 
from the observed distribution of planetary orbital periods. 
Subject headings: planetary systems — stars: individual (HD 37124) 



1. INTRODUCTION 

To date, over 50 exoplanetary systems with more than 
one planet have been discovered, including: the extraor- 
dinary detections of the firs t exo planets orbiting the pul- 
sar PSR B1257-hl2 (Wolszczan & Frail 1992; Wolszczan 
[1994); the imaged system orbiting HR 8799; those dis- 
cover eddimn^_tlwm^^ event OGLE-2006-BLC- 
109L (jCaudi et al.ll2008h : several systems discovered by 
transit, including four or fiveP^ multiply transiting sys- 
tems from the Kepler mission (Stcffcn 2010); and 43 sys- 
tems discovered by radial velocity (RV) searches (jWrightl 
l2009f l. The RV systems include the four-planet sys- 
tems u_AraJSantos ct 3] IIQOI jPeie^^ GJ 



581 jMavor et al.„2009i) and GJ 876 (jRivera et al.l 120051 
120101) and the five -planet system orbiting 55 Cancri 
([Fischer et al.l [2"008[ ). Of all these multplanet systems, 
only four are known to host three or more giani[^ plan- 
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et s with well-determined orbital paramaters: v And 
(Butl er et al . 1999), HIP 14810 (Wright ct al. 2009a), /i 
Ara ( Pepe et al.ll2007D and HD 37124' (|Vogt et al.l l2005). 

HD 37124 (HIP 26381) is a 85 metal- 
poor ('fFe/ H1=-0.44. iValenti fc Fisched [200l G4 dwarf 
(V=7.7). IVogt et all ([2OOO0 announced the a Jovian, 
P - 150 d planet orbiting HD 37124 from HIRES 
data taken at Keck Observatory as part of the Cali- 
fornia and Carnegie Planet Search. Further monitor- 
ing of the star rev ealed substantial long-term residuals. 
iButler et al.l ([2003| ) fit these residuals with an eccentric, 
1940 d planet, but noted that the solution was not unique 
(and' Gozdziewskil ([2003D showed that this fit was, in fact, 
unstable.) 

Afte r collecting two more years of data, IVogt et all 
l[2005h was able to report the detection of a third planet 
in the system, though with an ambiguity: while the b 
and c components had clearly defined periods, the d 
component could be fit nearly equally well with peri- 
ods of either 2300 d or 29.32 d, the latter hkely b eing 
an alias due to the lunar cycle [3 iWrightl ([20091 ) re- 
ported that recent Keck velocities had resolved the am- 
biguit y qualitatively in favor of the lon g er orb ital period. 
Gozd ziewski. Konacki. fc Macieiewskil ([2006D explored 
the many possible dynamical configurations consistent 
with the Vog t et al., (2005.) velocities, inclu ding many res- 
onant solutions. iGozdziewski. Breiter. fc B orczvk (200^ 
used the system to demonstrate a fast MENGO algo- 
rithm, but they did not explore the 2:1 resonance, as the 
data did not seem to favor it at the time. 

We present new Keck observations, and these data pro- 
vide a unique orbital solution for the outer planet. The 



1^ Time on the Keck telescopes dedicated to observing bright, 
planet search targets with HIRES is usually assigned during bright 
or gray time; the resulting scarcity of data points during new moon 
can interact with planetary signals to create spurious, aliased so- 
lutions. 
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outer planet period we find is more con s istent with the 
original period r eported bylButler et al.l ()2003l ) than the 
refined orbit of iVogt et alJ (|2005r FI (though we find a 
much lower eccentricity). Herein, we present the en- 
tire history of Keck velocities obtained for this star, and 
present self-consistent orbital solutions showing that the 
outer two planets are in or very near a 2:1 mean- motion 
resonance (MMR) . This is the 20th exoplanetary system 
to be found near an MMR, and only the tenth system 
with an apparent 2:1 commensurability. 

Period commensurabilities (PCs) represent impor- 
tant dynamical indicators in the Solar System and 
have been li nked with obse rvables and formation 
mechanisms (jGoldreichl I1965I ). The near-5:2 PC 
of Jupiter and Saturn, also known as "The Great 
Inequality", might be the remnant of a divergent 
resonant crossing that produced the current archi- 
tecture of the outer Solar System, the Late Heavy 
Bornb a rdment, and the T rojan Asteroids (iGomes et alJ 
20051: IMo rbidelli et a ll [2005 iTsiganis et alJ l2nn5r 



Tsiganis. V arvoglis. & DvorakI 20051 ). The populations 



of the asteroid belt and the Kuiper Belt, exemplified 
by t he PC and near-PC-populated Kirkwood Gaps 
(e.g. ITsiganis. Varvoglis, fc Hadiidemetrioul |2002|) the 
Plutinos (3:2 PCs with Pluto and Neptune) and 
the twotinos (2:1 PCs with Pluto and Neptune; e.g . 
IMurrav-Clav fc Chiand [2005t IChiang fc .ToTdm] l200l . 
have implications for the migratory history of Jupiter 
and Neptune and the prospect of, e.g. secular reso - 
nant sweeping (e.g., iNagasawa. Ida, fc Besshol I2008D . 
Near-PCs found in satellite and ring systems have 
had direct observational consequences; the Saturnian 
satellite Pandora was ~ 19° behind its pre dicted orbital 
longitude in a 1995 ring plane crossing (jFrench et al.l 
[2003) due to its 121:118 PC with neighboring satellite 
Prometheus. 

By extension, we may anticipate similar importance in 
the growing number of exoplanetary systems exhibiting 
PCs. In extrasolar systems, Mean Motion Resonances 
(MMRs) have been interpreted as the indication of 
convergent migration in multi-planet systems, (e.g. 



Thommes & LissaueHl2003 


:lKlev. Peitz. & Brvdenll2004l: 
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2005). Several subsequent 



Terouem fc Paoaloizou 2007; Picrcns & Nelson 2008, 



Radlcwska & Szuszkiewicz 2008, 2009; Li bert fc Tsiganis 
2009: .Rein fc Papaloizou 2009; Papaloizou fc Terquem 



20l3 iRein. PapaloYzou. fc Kiev! 120101; iZhang fc Zhou 
2010allbl ) exploring convergent migration for a variety of 



masses, separations and disk properties have found many 
regions of mass and orbital element phase space in which 
planets are easily captured through this mechanism. 

2. VELOCITIES AND ORBITAL SOLUTION 

Table [T] contains radial velocity m easurements for H P 
37124 from the HIRES spectrograph (Vog t et al.lll994l ) at 
Keck Observatory obtained by the Californi a Planet Sur- 
vey c onsortium using the iodine technique (jButler et al.l 
Il996f) . Note that the quoted errors are our internal (ran- 



IVoet et al.l ([2005') opted to refer to the new , 840 d signal as 
the c component, despite the prior 1940 d fit of Butler et al.': (2003), 
because that prior fit was so speculative, and because their new fit 
put the very existence of a 1940 d periodicity in some doubt. 



dom) errors, with no "jitter" included (jWrightl 120051 ). 

These velocities supersede our previously published ve- 
locities for this star, as we continue to refine our data 
reduction pipeline. Our ever-evolving radial velocity 
pipeline is descended in spirit and form from that de- 
scribed in Butler et al. (1996), but includes many small 
and large technical improvements, a thorough discussion 
of which is beyond the scope of this manusc r ipt. S ome 
details can be found in § 4.1 of IHoward et"aI1 (l20Tol) . § 3 
of lHoward et al.l (|2009| ). and in Batalha et al. 2011 (ApJ, 
accepted). 

On issue of instant relevance is that in August 2004 the 
HIRES CCD detector was upgraded to a CCD mosaic. 
The old Tektronix 2048 EB2 engineering-grade CCD dis- 
played a variable instrumental profile asymmetry due to 
a charge transfer inefficiency which manifested itself as 
small changes in a star's measured radial velocity as a 
function of exposure time (i.e. raw counts on the chip.) 
We apply an emperical, spectral-type dependent model 
to correct this effect for velocities measured prior to the 
detector change. The new CCD mosaic shows no evi- 
dence of this effect, but as a consequence of the switch 
there is a small velocity offset between data sets that span 
the two detector s ets similar to the detector -to-detector 
offsets discussed in lGregorv fc Fisclieil (|2010[ ). These off- 
sets could, in principle, be different for every target. 

Analysis of RV standards and known planetary systems 
show that such an offset is usually small - of order 5 
m/s - and very often consistent with zero. As a result, 
we report two independent data sets for this system in 
Table [1] one from each of the two detectors. We solve for 
the detector offset as an unconstrained free parameter. 
The times of observation are given in JD-2440000. 

We fitted the data using the publicly available multi- 
planet RV-fitting I DL p ackage RV_FIT_MP, described in 
I Wright fc HcwardI (pOOa ). In Table [2] we present our 3- 
planet Keplerian (kinematic) fitjH which yields r.m.s. 
residuals of 4.4 ms""'^, and we plot the fit and veloci- 
ties in Figure [T] We find a best-fit offset between CCDs 
of 4.8 ms""'^. The orbital parameters and their uncer- 
tainties were det ermined from 100 bootstrapped trials 
(as described in IMarcv et al.l I2005t IButler et all 120061; 
iWright et al.ll2007D . The orbital fits and dynamical anal- 
ysis herein are put forth under the assumption that the 
velocities are not detectably influenced by additional, 
unmodeled planets in the system. We have integrated 
these orbital parameters for 10 Myr using the methods 
described in 33] assuming coplanarity, and found them to 
yield a stable conflguration. 

The residuals to this flt have an RMS 4.03 m/s and 
show with no significant periodogram peak at any pe- 
riod. The tallest peak is at 3.81 days. We have run a 
Monte Carlo FAP analysis on these residuals of our best 
fit for this tallest peak, and find similarly good fits in 
50% of velocity-scrambled trials, consistent with noise. 
We thus conclude that our model is sufficient to explain 
the data and that there are no other statistically signifi- 
cant planetary signals detected. 

3. NEWTONIAN FITS AND STABILITY ANALYSIS 

This solution is of similar quality to the best fit Newtonian 
solution, and is dynamically stable. We consider it representative 
of the ensemble of good Newtonian solutions. 
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Possible 2:1 Resonance for HD 37124 
Radial velocity curves for the HD 37124 triple system. 
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TABLE 1 
Radial Velocities for HD 37124 



Time 


Velocity 


Uncertainty 


CC 


JD-2440000 


m 


m 




10420.04655 


54.81 


1.5 


1 


10546.73646 


28.74 


1.2 


1 


10837.76625 


6.32 


1.6 


1 


10838.94868 


6.38 


1.7 


1 


10861.80464 


17.89 


1.4 


1 


11069.03617 


-3.99 


1.4 


1 


11070.13190 


-1.66 


1.3 


1 


11071.11494 


1.28 


1.6 


1 


11072.12947 


-11.58 


1.5 


1 


11073.02962 


-8.72 


1.3 


1 


11172.89571 


39.40 


1.6 


1 


11226.78065 


-0.48 


1.4 


1 


11227.78167 


-2.05 


1.4 


1 


11228.74293 


-11.01 


1.3 


1 


11412.14161 


-33.76 


1.6 


1 


11543.98278 


-31.36 


1.4 


1 


11550.94262 


-44.57 


1.4 


1 


11551.94008 


-46.20 


1.5 


1 


11552.89162 


-47.97 


1.5 


1 


11580.76121 


-36.28 


1.8 


1 


11581.83559 


-36.17 


1.4 


1 


11582.78849 


-37.36 


1.4 


1 


11583.72387 


-35.69 


1.5 


1 


11884.04436 


-14.73 


1.6 


1 


11900.03518 


3.01 


1.4 


1 


11974.80019 


39.56 


1.4 


1 


12007.74522 


4.48 


1.5 


1 


12242.99064 


48.20 


1.5 


1 


12333.94545 


-19.13 


1.7 


1 


12334.78556 


-12.82 


1.7 


1 


12536.12848 


10.64 


1.8 


1 


12537.08597 


10.28 


1.7 


1 


12573.03767 


27.38 


1.6 


1 


12574.99934 


27.48 


1.7 


1 


12576.02212 


23.58 


1.6 


1 


12600.99996 


8.06 


1.7 


1 


12602.03213 


7.02 


1.7 


1 


12925.01639 


13.03 


1.7 


1 


13044.77359 


36.56 


1.6 


1 


13045.74638 


32.78 


1.5 


1 


13072.85948 


-2.85 


1.7 


1 


13240.13983 


-34.52 


1.5 


2 


13302.13193 


-4.28 


1.6 


2 


13302.97959 


-5.69 


1.4 


2 


13338.96446 


15.77 


1.2 


2 


13340.09520 


16.45 


1.5 


2 


13368.88930 


-6.51 


1.0 


2 


13369.78156 


-6.60 


1.0 


2 


13425.87197 


-39.98 


1.4 


2 


13426.82980 


-39.12 


1.3 


2 


13428.77030 


-36.85 


1.3 


2 


13483.72749 


19.49 


1.0 


2 


13723.90630 


-7.28 


1.6 


2 


13841.76427 


34.09 


1.4 


2 


14544.83023 


25.40 


1.6 


2 


14545.78169 


26.99 


1.4 


2 


14546.78977 


24.28 


1.3 


2 


14718.08322 


43.60 


1.6 


2 


14806.91704 


-7.74 


1.5 


2 


14810.89031 


-9.50 


1.7 


2 


14838.94681 


0.39 


1.8 


2 


14864.95362 


29.45 


1.8 


2 


14929.76349 


-5.49 


1.7 


2 


15135.00085 


-21.29 


1.6 


2 


15172.99171 


12.60 


1.6 


2 


15229.78574 


-12.32 


1.6 


2 



TABLE 2 

Best-Fit Kinematic Orbital Elements for Exoplanets in the HD 
37124 System 



Parameter b c d 



P(d) 


154.378 ± 0.089 


885.5 ± 5.1 


1862 ± 38 


Tp (JD-2440000) 


10305 ± 11 


9534 ± 11 


8558 ± 11 


e 


0.054 ± 0.028 


0.125 ± 0.055 


0.16 ± 0.14 


CO (degrees) 


ISO^' 


53 ± 17 


Qa 


K (ms-l) 


28.50 ± 0.78 


15.4 ± 1.2 


12.8 ± 1.3 


Afsinj(Mjup) 
a (AU) 


0.675 ± 0.017 


0.652 ± 0.052 


0.696 ± 0.059 


0.53364 ± 0.00020 


1.7100 ± 0.0065 


2.807 ± 0.038 


R.M.S. 




4.03 




jitter 




4 m/s 








0.8 




^obs 




66 





^ Orbit is consis tent with circular, so errors in lu are large; see (see 
IButler et al.ll200^ . for a fuller explanation.) 



3.1. MCMC analyses 

We studied the dynamical stability of HD 37124 by 
combining the radial velocity data with Markov Chain 
Monte Carlo (MCMC) analyses to obtain ensembles of 
masses, semimajor axes, eccentricities, and orbital an- 
gles consistent with the RV data. These ensembles were 
generated without regard to dynamical stability consid- 
erations. We then imposed line-of-sight and relative in- 
clination distributions on these sets of parameters. By 
incorporating the unknown inclination parameters with 
observation-derived parameters, we sampled the entire 
phase space of orbital parameters. We subsequently ran 
N-body simulations on each element in these ensembles 
in order to assess each system's stability and re s onant 
evolu t ion. Our treatment follow s that of iFordI ()2005l 
l20nfih : l\^s fc FordI (|2?M l20Tl . 

In particular, we calculated 5 Markov chains, each con- 
taining over 10^ states. Each state includes the orbital 
period (P), velocity amplitude (K), eccentricity (e), ar- 
gument of pericenter measured from the plane of the sky 
(w), and mean anomaly at a given epoch (u) for planets b, 
c and d. The MCMC uses a standard Gaussian random 
walk proposal distribution and the Metropolis-Hastings 
algorithm for accepting or rejecting each proposal for all 
model parameters except cos(iLOs) and Q. Since the ra- 
dial velocity signature is only weakly dependant on these 
values, cos(iLOs) Sind f2 were drawn randomly for each 
state. This can still be considered a Markov chain, as the 
procedure satisfies the Markov condition, i.e. that a trial 
state not depend on states other than the current state, 
as well as the other conditions (time-homogeneous, irre- 
ducible, aperiodic) to prove that the Markov chain will 
(eventually) converge to the posterior distribution. 

We imposed an isotropic distribution of line-of-sight 
inclinations (ihos) and a uniform sample of longitude of 
ascending nodes (fi) on our MCMC-derived initial con- 
ditions. The planet masses, m, and semimajor axes, 
a, were obtained from each set of (P, if, e, oj, i, $7, it) 
values fr om relations derived with a Jacobi coordinate 
system (|Lee fc Peald 120031 ). The approximate range 
of minimum masses obtained, in Jupiter masses, were: 
0.60 < mbsinifc < 0.72, 0.40 < rricSmic < 0.75, and 
0.55 < mrfsinirf < 0.90. 

We treated the both the offset between the chips and 
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the jitter as free parameters The 5 percentile, me- 
dian, and 95*^ percentile offsets between the chips in our 
ensemble were 3.16, 3.78, 4.62, and the median jitter we 
find to be 4 m/s. 

3.2. Coplanar, Prograde Systems 

We integrated 850 sets of initial conditions in the copla- 
nar case with all three planets in prograde or bits by us- 
ing the Burlish-Stoer integrator of Mercury (jChamberd 
Il999f l for 10^ yr with an output interval of lO'' yr. 
We also incorporated the effects of general relativity in 
the code, which could have p rofound consequence s for 
multi-planet system stability (jVeras fc FordI l20To() . al- 
though the effect is likely to be negligible in this system. 
We classified systems as "unstable" if, for any planet, 
Iflmax - flminl/ao > T, where amax, Omin and ao repre- 
sent the maximum, minimum and initial values of the 
semimajor axis, and t = 0.9. 

One may visualize a representative architecture of HD 
37124 by comparing the semimajor axis and eccentricity 
ranges of all three planets. Figure [2] plots the observed 
eccentricity vs. derived semimajor axis for all planets in 
the prograde coplanar state. Black dots indicate unstable 
systems while green squares and red crosses indicate sta- 
ble systems, and red crosses indicate systems which are 
in a 2:1 mean motion resonance (MMR) between planets 
c and d, according to our definition below. The figure 
indicates i) a closely packed system, with the inner and 
outer planets separated by no more than six times the 
innermost planet's semimajor axis, ii) a relatively circu- 
lar innermost planet (with e;, < 0.1 in most cases) that 
is likely too far from the parent star to be classified as 
a "Hot Jupiter", iii) the greater the number of orbital 
periods sampled by RV, the greater the constraint on 
the planet's likely semimajor axes and eccentricities, iv) 
most (664/850 = 78%) current orbital fits predict an un- 
stable system, v) the majority of initial conditions which 
produce stable orbits contain an outer planet with a low 
(< 0.2) eccentricity, and a middle planet with a semi- 
major axis > 1.695 AU and eccentricity less than about 
O.2., vi) systems containing a 2:1 resonance occur only 
when 2.7 AU < 03 < 2.8 AU. We emphasize that this 
approximate semimajor axis range appears to be neces- 
sary but not sufficient for resonance to occur. The figure 
demonstrates that other MCMC fits with outer planet 
periods in the resonant range are either unstable, or sta- 
ble but non- resonant. The architecture of these systems 
(as defined by, e.g., the mean longitude and longitude 
of pericenter) do not allow them to settle into resonance 
even though the outer planet period might favor reso- 
nance. 

Because of finite sampling, our definition of "reso- 
nance" in this analysis comes from consideration of the 
RMS deviation of each resonant angle about each of 
(0°, 90°, 180°, 270°), which includes common libration 
centers. We flag systems as "resonant" if at least one 
of these angles has RMS under 90° for 10 Myr, the en- 

We adopted a single value of jitter for all observations; in prin- 
ciple the two CCDs may display differing amounts of "instrumen- 
tal jitter", sucli as tliat due to insufficient modeling of the charge 
transfer inefficiencies. The RMS residuals to our fit for the two 
detectors were 3.67 and 4.11 m/s, suggesting that our assumption 
of a single jitter value is valid. 



tire duration of our simulations. Below, we refer to this 
value as a "libration RMS" . 

HD 37124 presents a clear initial choice of angles to 
test for libration. As indicated by Fig. [21 the semimajor 
axis ratio of planets c and d is suggestive of a 2:1 MMR. 
Therefore, we sampled the following angles for libration: 

01 = 2Arf - Xc- vjc (1) 

02 = 2A<j - \c~ r^d (2) 

and found that 01 librates in 28/850 — 3.3% of cases, 
while 02 librates in 9/850 = 1.1% of cases. Further, the 
systems for which 02 is resonant are a subset of those for 
which 01 is resonant. 

If we tighten our definition of resonance to include 
only those systems with RMS resonant angles under 
70°, then no 02 arguments are resonant. Under this 
stricter definition, the 0i arguments are only resonant 
in 14/850 = 1.6% of cases, and if we further tighten the 
libration criterion to an RMS of or 50°, then this num- 
ber decreases to 4/850 — 0.5%. The lowest libration 
RMS detected is 23.0°. AU RMS's under 75° were for a 
libration center of 0°. Figure [3] illustrates three examples 
of "resonant" systems from this, each with a different li- 
bration RMS. 

We additionally sampled all 3 pairs of apsidal angles 
(the difference between two longitudes of pericenter) in 
the coplanar prograde state, and found only two in- 
stances of libration, both at high (> 70°) libration RMS's 
and around the "asymmetric" centers 90° and 270° for 
the inner and outer planet apsidal angle. Inspection re- 
veals, however, that these instances of libration are more 
indicative of long period (> 10 Myr) circulation. 

Additionally, the semimajor axis ratio of planets b and 
c could indicate the presence of a 6:1 MMR. Therefore, 
we sampled all angles of the form 6Xd — Ac — twc — swd, 
where t + s — 5. None of the coplanar prograde sys- 
tems exhibited libration of any of the 6:1 angles between 
planets "b" and "c" over 10 Myr. However, preliminary 
sampling of these angles over intervals of 2 Myr does oc- 
casionally exhibit libration RMS's close to 90°. Because 
the period ratios between planets c and d may skirt the 
7:3 PC, we also tested the 7Xd — 3Ac — twc — swd angles, 
where t -\- s = A, but found no instances of libration. 

3.3. Mutually Inclined Systems 

Having analyzed the coplanar prograde bin, we can 
now consider the case where the planets have nonzero 
mutual inclinations. We used rejection sampling to ob- 
tain triplets oil LOS values such that the system is placed 
into one of 144 bins according to the relative inclination 
between planets b and c {irei,b,c) and planets c and d 
{irei,c,d)- In no two bins were the same ensemble of ini- 
tial conditions used. We binned relative inclinations in 
intervals of 15°, and used stratified sampling in order to 
obtain a uniform number of samples per bin. Wc initially 
sampled 100 initial system states per bin. For those bins 
where we found more than one system to be stable, we 
added 200 additional ensembles of initial conditions 

By considering the fraction of stable systems in the 
non-coplanar cases, we can obtain a broader dynamical 
portrait of this system. Fig. [4] illustrates the fraction of 
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Fig. 2. — Representative eccentricities and semimajor axes of the planets in HD 37124. The three planets are partitioned by panels, 
each with a different horizontal scale. These ensembles of parameters are derived from RV observations using Markov Chain Monte Carlo 
(MCMC) techniques, and represent the initial conditions for a subset of our numerical simulations (here, the coplanar prograde simulations). 
Note that the semimajor axis ratio of the middle and outer planets roughly correspond to a 2:1 period commensurability, and the inner 
and middle planets to a 6:1 PC. Note the changing scale on the x-axes in the three panels: the innermost planet is very well constrained in 
a, the outermost planet much less so. The Black dots indicate unstable systems, green squares represent non-resonant stable systems, and 
red x's are stable resonant systems. Note that system stability is strongly dependent on the eccentricity of the middle and outer planets, 
and the outer planet of resonant systems tend to harbor the smallest initial semimajor axes of the ensemble of outer planet ICs. 



stable systems in each bin overall (top panel) and with 
respect to all systems for which the initial < 0.2 (bot- 
tom panel) . This cutoff was motivated by the rightmost 
panel in Fig. [3] and could suggest a constraint on the 
orbital properties of the system in order to ensure that 
it remains stable. Fig. 2] shows that the system must 
be roughly coplanar, with relative inclinations less than 
~ 30° — 45°, in order to be stable. This constraint al- 
lows various pairs of planets to harbor retrograde orbits. 
We also performed limited resonant testing for systems 
in these bins. The fraction of total systems which exhibit 



libration of 01 and 02 under 90° for 10 Myr is given by 
Fig. [5l 2:1 resonant systems occur, therefore, generally 
at the few percent level, and most likely when all planets 
are coplanar with prograde orbits. 

4. DISCUSSION 

As lRivera et al.l (j2010) showed for the GJ 876 system, 
even the most well-established and deepest mean-motion 
resonances can prove illusory if additional planets are 
found in the system (although in that case it appears 
that the resonance still present, albiet considerably shal- 
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Fig. 3. — Three examples of systems that we find to be resonant, according to our definition requiring an libration of under 90° for 10 
Myr. Plotted is the time evolution of the resonant angle 2X^ — Ac — ■cUc for a system with a computed libration RMS of (upper panel) 23.0° 
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lower and more complex than previously thought). Even 
for truly resonant systems, a demonstration of resonance 
can be difficult. For instance, triple-planet systems may 
feature two planets with a mostly-librating resonant ar- 
gument that occasionally circulates due to interactions 
with the third pl anet. Near separatrix behavior (as in 
the case of v And; lMalhotrall2002l: iFord. Lvstad, fc Rasiol 
[20051) can also make libration and circulation essentially 
indistinguishable . 

We note that near-resonant behavior can itself dynam- 
ically interesting: the 5:2 near-resonance of Jupiter and 
Saturn (the Great Inequality) has major consequences 
for the dynamics of the Solar System. Given the above- 
mentioned difficulties in proving that a resonant argu- 
ment for a given system of planets satisfies some pre- 
cisely specified definition of libration given the typical 
uncertainties in radial velocity measurements, we sug- 
gest that studies of resonant interactions would benefit 
from identifying systems that appear to be in or near res- 
onance {apparent period commensuribilities) . With that 
in mind, we note that in addition to HD 37124, there 
are 21 other systems in the Exoplanet Orbit Database 
(Wright et al. 2010, PASP submitted) of peer-reviewecO 
literature with well-established apparent PCs, which we 
present in Table |31 This list includes all pairs of plan- 
ets for which the period ratio r is less than 6 and within 
0.05 of an integer or half integer (neglecting uncertainties 
in periods), and other exoplanetary pairs whose likely 
MMRs are discussed in the literature. 

The fraction of known multiplanet systems exhibiting 
at least one apparent PC is high. Of the 43 well de- 
termined multiplanet systems discovered by radial ve- 
locities around normal stars, 15 appear in Table [3l or 
35%, including 9 of the 30 apparent double-planet sys- 
tems 30% □ To determine if this is more than would be 
expected simply by chance, we have performed two tests. 

In the first test, we randomly drew pairs of periods 
from the 340 RV-disc overed planets in the Exoplanet Or- 
bit Database (EOD) ([Wright et al.ll2010l PASP submit- 
ted), and rejected those pairs with period ratios r < 1.3 
(corresponding to the smallest r among real multiplanet 
systems) . We counted the fraction of remaining systems 
with r within 0.05 of an integer or half integer < 5 (cor- 
responding to the largest apparent PC in Table [3|). We 
found that only 4% of our random pairs satisfy our ap- 
parent PC criterion, far smaller than the 30% of double- 
planet systems actually found in apparent PCs. 

In the second test, we included the effects of triple and 
higher-multiple systems by randomly assigning periods 
from the EOD to all planets in real multiple systems 

^® We have included the Kepler multiplanet systems, which had 
not yet been accepted for publication at the time of this writing. 

We have excluded in this statistic the Kepler systems, the 
pulsar system, the microlensing system, planets from direct imag- 
ing, and the Solar System. We acknowledge that a more rigor- 
ous statistic would be valuable, but note that it would need to 
address some strong detectability and selection effects regarding 
planets in multiple planet systems, and to assess these detection 
thresholds across multiple, heterogeneous surveys. To give just 
one example, we note that in addition to a radial velocity sur- 
vey's decreasing sensitivity to planets in longer periods, it can 
be difficult to detect an interior planet in a 2:1 resonance due 
to approximate degenerac y with eccentricity in a s ingle planet 
model (e.g., A nglada-Escu de. Lopez-Morales, fc Cha mbers 2010; 
IMoorhead &: PordI IZOIOD . Such an analysis is beyond the scope 
of this manuscript. 



(again subject to the constraint that no pair of planets 
in the system have r < 1.3). We found 16% of these 
artificial systems passed our apparent PC criterion, re- 
flecting the higher number of planet pairs available to 
test per star compared to our first test. Despite this in- 
flation, the actual value of 33% among all multiplanet 
systems is still significantly higher. [3 These results 
underscore that the orbital periods of the population of 
planets known to be in multiplanet syste ms is inconsis- 
tent with the apparently-singleton sample ([Wright et al.l 
2009b). 

This apparently high percentage of known multi- 
planet systems with an apparent PC might favor par- 
ticular formation mechanisms. Planet-planet scatter- 
ing, planetesimal disk migration, and gas disk migra- 
tion have all been shown to produce systems with at 
least one pair of planets that n ot only are commensu - 
rate in period, but also resonant. IRavmond et al.l ([2008f ) 
found that planet-planet scattering produced MMRs in 
roughly 5% of the systems tha t they simulated, and 
IRavmond. Armitage. fc Gorelickl ([2010D discovered that 
between 50% - 80% of systems undergoing planetesimal 
disk migration yielded resonant capture. Convergent 
gas disk migration, the thrust of the numerous papers 
cited in [JT] can occur with near 100% efficiency for cer- 
tai n initial planetary and disk p aram eters. As observed 
bv [Th ommes fc Lissaueil ([2003h and [Libert fc TsiganisI 
(|2009( ), the inclination may be excited as well as the ec- 
centricity in many resonant cases. If a resonance exists in 
HD 37124, it could have been produced by any of these 
methods. If the RMS libration of such a resonance is rep- 
resentative of the bottom panel of Fig. [3| and has a value 
that approaches 90 degrees, then planet-planet scatter- 
ing is a likely origin of this resonance. Alternatively, disk 
or gas migration would likely produce a system that is 
"deeper" in resonance, with a smaller variation in reso- 
nant angle, similar to the top panel of Fig. [3] Resonant 
librating angles need not involve the eccentricities and 
pericenters, but instead the inclinations and longitudes 
of ascending nodes, similar to the 4:2 Mima s- Tcthys res- 
onance in the Saturnian system ([Champenois fc Vienn^ 
il999bi ial). 

5. CONCLUSIONS 

We have resolved the period ambiguity of HD 37124 
d from [Vogt eFall ([2005[ ) and find that HD 37124 c and 
d are in an apparent 2:1 period commensurability. Our 
numerical integrations show that both resonant and non- 
resonant configurations are consistent with the radial ve- 
locity data, and that stability requires a nearly circular 
orbit (e < 0.3) for the d component. Our stability anal- 
ysis shows that the system must be nearly coplanar, and 
that the three planets have identical minimum masses 
within the errors (of 3-10%). 

We show that the roughly one in three well- 
characterized multiplanet systems shows an apparent pe- 
riod commensurability, which is more than a naive esti- 
mate based on randomly drawing periods from the known 
exoplanet population would suggest. This offers evidence 

We suspect that the reason the observed value is not similarly 
inflated with respect to double planet systems is that our random- 
ization did not include the requirement of dynamical stability, as 
real systems implicitly do. 
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TABLE 3 

Apparent period commensuribilities in well-characterized multiplanet systems 



System 



Components Period Ratio 



References 



GJ 876 


G c b 


4:2:1 


HD 82943 


b,c 


2:1 


HD 37124 


c',d 


2:1 


HD 128311 


c,b 


2:1 


HD 73526 


c,b 


2:1 




b,e 


2:1 


KOI 152t 


2,3 


2:1 


KOI 877t 


2,1 


2:1 


24 Sex 


c, b 


2:1 


KepIer-9 


c,b 


2:1 


PSR B1257+12 


B,C 


3:2 


HD 45364 


c,b 


3:2 


HD 200964 


c, b 


4:3 


55 Cnc 


c, b 


3:1 


HD 10180 


d, e 


3:1 


HD 60532 


c, b 


3:1 


HD 108874 


c, 6 


4:1 


Solar 




5:2 


HD 10180 


e,f 


5:2 


KOI 896t 


1,2 


5:2 


HD 202206 


c, 6 


5:1 



Marcv et al. (20011 : Rivera et al. (2010) 

Mayor et al. (2004): Ji et al. (2003TTee et al. (2006) 

Vogt et al. (2005), this work 

Sandor & Kley (2006a); Vogt et al. (2005), others 
Tinney et al. (2006); Sandor & Kley (2006b) 

^e£e_et_aL (_2007 ) ; ^ozdxiewskij_^aciejewBkij_^^igas^ (_2007 ) 
Steffen (2010) 
Steffe n (2010) 
Johnson et al. (2011) 
jlolman et al. (2^To) 

Wolszczan & Frail (1992); Malhotra et al. (1992) ; Konacki. Macieiewski. 
Rein, Papaloizou, fc Kiev (2010) ; CorTeiTet al. (2009) 
Johnson et al. (20n) 

Fische^_eta^ . . , ^^^^^^^ 

Lovis et al. 201CirA&T accepted 
Desort et al. (2008 . 2009 ); Laskar 
^og^_et_^l^ 



Wolszczan (1999 




Correia (2009) 

Maciejewski 



Lovis et al. 2010, A&A accepted 
Steffen (2010) 

^orrei^_et_aL (2005) ; ^ozdaiewskij_j<on^ddj_^^ac^^ (2006) 



t This is a candidate exoplanet system based on Kepler photometry, but the planets are at present considered unconfirmed. For example, KOI 
877 could, in principle, be a blend of two stars, each hosting one transiting planet in a coincidental apparent PC. 



for some particular proposed scattering and migration 
mechanisms, and we suggest that the statistics of multi- 
planet systems may now be sufficiently robust to provide 
a test and comparison of models of exoplanetary dynam- 
ical evolution. 
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